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Motivated by a recent experiment on Rb2Cu3SnFi2, where spin-1/2 Cu^^ moments reside on the layers of 
Kagome-like lattices, we investigate quantum ground states of the antiferromagnetic Heisenberg model on a 
series of deformed Kagome lattices. The deformation is characterized by a weaker exchange coupling a J on 
certain lattice links appropriate for Rb2Cu3SnFi2 with a — 1 corresponding to the ideal Kagome lattice. In par- 
ticular, we study possible valence bond solid phases using the perturbation theory around isolated dimer limits, 
dimer series expansion, and self-consistent bond operator mean field theory. It is shown that the valence bond 
solid phase with a 36-site unit cell of the ideal Kagome lattice is quite sensitive to a small lattice distortion as the 
kind discovered in Rb2Cu3SnFi2. As a result, we find that a more likely quantum ground state in Rb2Cu3SnFi2 
is the valence bond solid phase with a 12-site unit cell, where six dimers form a pinwheel structure, leading to 
strong modification of the elementary triplet and singlet excitation spectra in the deformed Kagome lattices. 
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I. INTRODUCTION 

Lattice distortions are abundant in frustrated magnets as 
it is an efficient way to relieve the frustration and lift the 
macroscopic degeneracy of the classical ground states. This 
has been, for example, an obstacle to synthesize and study 
materials that can be described by the quantum Heisenberg 
model with spin-1/2 moments on "ideal" lattice structures 
with isotropic spin exchange couplings. The interest on the 
"ideal" models stems from various theoretical suggestions for 
novel emergent quantum paramagnetic phases such as quan- 
tum spin liquid''^'^'''*'^'* and valence bond solid (VBS)^''*'^'*^^ 
phases. On the other hand, those deformations may lead to 
different kind of emergent quantum ground states and un- 
derstanding precisely how these phases are connected to the 
quantum ground states in the "ideal" limit would help us 
understand more global picture of the quantum frustrated 
magnets."''^ 

One of the most studied frustrated magnets is the quan- 
tum spin-1/2 antiferromagnetic Heisenberg model on the 
Kagome lattice. There exist several materials that are can- 
didates for the realization of this system. '^''''''^i'^''^''^i'VO 
It has been, however, found that many materials undergo 
lattice deformations and quite often this leads to magnet- 
ically ordered phases at sufficiently low temperatures.'^''"' 
So far, there are a small number of examples of two- 
dimensional Kagome or Kagome-like materials with spin- 
1/2 moments that do not show any magnetic ordering 
down to very low temperatures. '^''^''^''^''^'^° The Herbert- 
smithite, ZnCu3(OH)6Cl2 ,'^i'^'"^ is known for the "ideal" 
Kagome lattice structure, but it suffers from inter-site mix- 
ing between Cu-^+ and Zn^+, and makes the interpreta- 
tion of the susceptibility data difficult. In the Volborthite, 
Cu3V207(OH)2-2H20,''^ the Kagome structure is deformed 
in an orthorhombic fashion, leading to two inequivalent spin 
exchange couplings between Cu^+ moments. In both cases, 
the specific heat at low temperatures reveals the absence of 
(or a very small) excitation/spin gap (down to 50mK in Her- 
bertsmithite, for example), leading to various suggestions for 
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FIG. 1: (Color online) (a) The deformed Kagome lattice with four 
different neighboring exchange interactions, which may be realized 
in Rb2Cu3SnFi2. Here all the exchange coupling constants are pos- 
itive and satisfy Ji > J2 > J3 S> J4. In this work, we consider the 
case of Ji = J2 = J3 = J > J4. (b) The 14-depleted Kagome lat- 
tice that is topologically equivalent to the deformed Kagome lattice 
when Ji — J2 = J3 ~ J but J4 — 0. 



exotic quantum ground states.^'*-'' More recently, it is found 
that Cu^+ moments in Rb2Cu3SnFi2 form the layers of a de- 
formed Kagome lattice and do not show any magnetic order 
at low temperatures.^" The susceptibility data clearly shows 
a spin excitation gap (in contrast to the previous two materi- 
als) which is of the order of 20K while the average strength of 
the exchange coupling is estimated to be around 200K, which 
leads to the expectation that the ground state may be a VBS 
state. 

In this paper, we investigate the quantum ground states 
of the spin-1/2 antiferromagnetic Heisenberg model on a se- 
ries of deformed Kagome lattices that may be relevant to 
Rb2Cu3SnFi2. Each deformed Kagome layer can be de- 
scribed by four possibly different nearest-neighbor exchange 
couplings, labeled by Ji, J2, J3, and J4 in decreasing order in 
magnitude as inferred from different Cu^+ - - Cu^+ bond 
angles (See Fig.IlJa)). In Rb2Cu3SnFi2, Ji, J2, and J3 are 
similar in magnitude while the magnitude of J4 is about a half 



of the others.— Hence it is reasonable to start from a model 
where Ji = J2 = J3 = J 3> J4. Here we study the evolution 
of the quantum ground states as a function of a = J4/ J < 1, 
where a = 1 corresponds to the ideal Kagome lattice. 

Using perturbative arguments starting from isolated dimer 
limits, we find that there exist two competing VBS phases on 
these lattices; one with a 12-site unit cell (VBS-12, shown in 
FigEJa)) and another with a 36-site unit cell (VBS-36, shown 
in Fig|3lb)). The VBS-36 phase is smoothly connected to the 
VBS phase with a large unit cell on the ideal Kagome lattice, 
which has been suggested as a strong contender for the quan- 
tum ground state of the spin- 1/2 antiferromagnetic Heisen- 
berg model. ^•'^•'^ The perturbation theory clearly shows that the 
VBS-12 is the lower energy ground state when a ^ 0, namely 
on the J4-depleted lattice. On the other hand, the VBS-36 
becomes the ground state when a ~ 1, consistent with the 
previous studies on the ideal Kagome lattice. We study the 
relative stability of these phases for an arbitrary deformation 
parameter a using the perturbation theory, dimer series expan- 
sion, and bond operator mean field theory. Thus our work can 
also test the stability of the VBS phase with a 36-site unit cell 
on the ideal Kagome lattice against the deformation from the 
ideal Kagome lattice structure. 

It is found that the VBS-36 phase on the ideal Kagome lat- 
tice is quite sensitive to the deformation and only 3% reduc- 
tion in the magnitude of J4 is enough to induce the instability 
of the VBS-36 state of the ideal Kagome lattice. Thus we 
propose that the ground state for a < 0.97 is the VBS-12 
state. Note that the VBS-12 state is clearly the ground state of 
Ji-only model or when J2, 0/3 ^ Ji. It is remarkable that the 
VBS-12 state is stable even when the strength of J2 and J3 are 
equal to Ji. Thus we conclude that the VBS phase discovered 
in Rb2Cu3SnFi2 is most likely the VBS-12 state. In order to 
provide a future experimental test for the VBS-12 state, we 
compute the triplet dispersion spectra as shown in Fig|6] 

The rest of the paper is organized as follows. In SecUll 
we use perturbative arguments about isolated dimer limits to 
show that there exist two competing VBS phases - the VBS-12 
and VBS-36 - on the J4-depleted lattice. In Sec|III] the dimer 
series expansion and the bond operator mean field theory are 
developed to investigate the relative stability of the two candi- 
date VBS phases for an arbitrary strength of J4/ J. The triplet 
dispersion spectra is also discussed in detail for the VBS-12 
state. Finally, in SeclIVI we discuss the implications of our 
results to theory and experiments. 



II. VALENCE BOND SOLID GROUND STATES ON THE 
J4 -DEPLETED KAGOME LATTICE 

A. Second Order Perturbation and the VBS-12 State 

We can obtain valuable information about the spatial or- 
dering pattern of singlet dimers by applying the perturba- 
tion theory around an isolated dimer limit. ^' To perform the 
calculation, we divide the Hamiltonian into two parts, i.e. 
H = iJo + . Here Hq describes the Heisenberg spin ex- 
change interaction between two nearest-neighbor spins con- 
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FIG. 2: (Color online) Local configurations of neighboring dimers, 
which contribute to the lowest order perturbation theory. Here an 
eUipse indicates a dimer singlet. A solid line shows the intradimer 
coupling included in Hq while a dotted line indicates the interdimer 
coupling contained in the perturbation, V . (a) Neighboring dimers 
are lying in parallel. There are energy gains Ae''^' = — ^ JA^ in 

the second order and Aei''' = ~-^J^^ in the third order, (b) 
Neighboring dimers are lying in a perpendicular geometry. There 
is no energy gain in the perturbation theory (explained in the text). 
Here the triangle carrying a dimer is called 2i filled triangle, (c) An 
empty triangle. There is no energy gain in the third order perturba- 
tion theory, (d) A perfect hexagon. There is an additional energy 
gain Ae|,^' = —-^JX^ in the third order perturbation theory. 



stituting a singlet dimer The inter-dimer interactions are con- 
tained in V which is treated as a perturbation. The ground 
state of the unperturbed Hamiltonian Hq, i.e. I^Po) can be 
taken as the direct product of a set of isolated dimer singlets 
lying on the J4-depleted Kagome lattice. 

Let us start from two possible local configurations of two 
neighboring dimer singlets as shown in Fig.|2a) and (b). Note 
that a pair of neighboring dimers can be aligned in parallel, 
connected by a single unoccupied link (Fig. [2a)) or they can 
be aligned in a perpendicular direction if one end of a dimer 
is simultaneously connected to both spins of its neighboring 
dimer (Fig. |2jb)). The important point is that the ground 
state energy can be lowered through second order perturba- 
tion processes only when the two dimers are lying in parallel. 
When two dimers are lying in a perpendicular direction, V 
|4'o)=0 because of the odd parity of the singlet state with re- 
spect to the reflection and there is no energy shift due to the 
second order processes. The energy gain from each configu- 
ration where two neighboring dimers are lying in parallel is 

In order to determine the global structure of the VBS 
phases, however, we need to take into account all the geomet- 
ric constraints for possible VBS patterns. Note that the basic 
building blocks of the J4-depleted Kagome lattice are trian- 
gles, hexagons and bridge links which connect each triangle 
to its neighboring hexagons (shown in Fig. [Ub)). One can 
show that the number of triangles, N/\, is given by N/\ = 
where N is the number of spins (or lattice sites) on the lat- 
tice. Similarly, the number of hexagons, A^Hcxagon, and the 
number of bridge links, iVeridgo, satisfy A^Hoxagon = f| and 
A^Bridgc = -y, respectively. The total number of links, iLink, 
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on the lattice is then given by 



^Link = 37Va + 67VHexagon + A^Bridgo = 

Since the number of dimers, A/^Dimci, is equal to ^ and each 
dimer is lying on a single link, the total number of the unoc- 
cupied links is iVLink — Nuima — — = N. The un- 
occupied links connecting neighboring dimers in various local 
dimer configurations are shown as the dotted lines in Fig.|2] 

According to the result of the second order perturbation the- 
ory for two neighboring dimers, when one of the links of a 
triangle is occupied by a dimer singlet (we call such a triangle 
as a filled triangle), the other two links would generate no en- 
ergy gain from the second order processes because these two 
unoccupied links are simultaneously connected to one end of 
a nearby dimer, generating a 'perpendicular' arrangement be- 
tween the filled link of a triangle and a nearby dimer Hence 
if we denote the number of the filled triangles as A^fiUcdA, 
the number of the unoccupied links that connect two nearby 
dimers in parallel configuration is A^Link— ^Dimoi — 2A^finodA- 
Therefore, the total energy gain per spin from the second order 
processes is given by 



{■ 



vfillcdA 



N 



(1) 



Note that the condition AffinodA=0 uniquely determines the 
ground state VBS configuration which is shown in Fig. Oa). 
Here the unit cell consists of six dimer singlets which lie in a 
pinwheel structure and we call this VBS state as the VBS-12 
state. In the VBS-12, every triangle is an empty triangle, that 
is, no link of a triangle is occupied by a dimer singlet. This 
is the dimer configuration one would obtain when only Ji is 
present or J2, J3 ^ Ji (See Fig.lTJa)). It is interesting to note 
that the VBS-12 state is stable all the way to the limit where 
J2 and J3 become comparable to Ji . 



B. Third Order Perturbation and the VBS-36 State 

Now we extend the perturbation theory to the third order in 
A. Let us first consider two neighboring dimers lying in par- 
allel as shown in Fig|3a). In addition to the energy gain from 
the second order perturbation, Ae^^^ this local configuration 
allows a third order process which results in the energy gain, 
Aei^'' = —^gJ^^- On the other hand, the dimer configura- 
tion in Fig|2|b) does not allow any third order process. In fact, 
this configuration is locally an eigenstate of the Hamiltonian 
and does not produce any virtual state in any order of the per- 
turbation theory. Similarly, the unoccupied links in Figllfc) 
making an empty triangle, would not generate any energy gain 
in the third order, either. Therefore when we calculate the 
third order energy gain coming from various local configura- 
tions, we have to remember that the unoccupied links on every 
filled and empty triangles would not produce any third order 
energy shift. 




FIG. 3: (Color online) Two candidate valence bond solid configu- 
rations on the deformed Kagome lattice. The figure on the left-hand 
side is the configuration on the original deformed Kagome lattice and 
the picture on the righthand side represents a topologically equivalent 
lattice of the J4-depleted lattice. The light dotted (yellow) lines on 
the lefthand-side pictures show the J4 links, (a) VBS-12: A valence 
bond solid with a 12-site unit cell, (b) VBS-36: A valence bond solid 
with a 36-site unit cell. Here the thick dotted (pink) lines indicate the 
unit cells. 



FiglHd) shows another configuration in which a new third 
order process can be generated. This structure, which is com- 
monly known as a perfect hexagon, consists of three neigh- 
boring dimers lying on a hexagon. Through the resonating 
third order processes around a perfect hexagon, the ground 



state energy can be lowered by Ae 
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Based on these results, the total energy gain per spin by the 
third order processes can be summarized as follows. 



A£;(3) =A6(3){-^i^-'^ 
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where A'^poi fcct represents the number of perfect hexagons and 
we use the fact that the number of triangles, N/s = 

satisfies Na = A^cmptyA + A^fillcdA- 

Therefore if we can maximize the number of the filled trian- 
gles and perfect hexagons, we can get the largest energy gain 
from the third order processes. The VBS-36 state depicted 
in Figl3jb) satisfies this requirement, that is, A'finodA — -^a 
and iVpcrfcct = I^^Hoxagon- Notc that the three neighboring 
hexagons lying in a triangular structure cannot be 'perfect' at 
the same time for any kind of dimer coverings. Interestingly, 



this VBS-36 phase is adiabatically connected to the 36-site- 
unit-cell VBS ground state of the ideal Kagome lattice**'^ in 
the sense that as we increase the magnitude of J4/J from 
zero to one, the VBS-36 phase evolves smoothly to the 36- 
site-unit-cell VBS ground state of the ideal Kagome lattice. 

In the case of the J4-depleted lattice, the VBS ground state 
is already uniquely determined by the second order perturba- 
tion theory. That is, the VBS-12 phase is the ground state and 
the VBS-36 phase is a closely competing phase. However, it is 
not clear which phase is more stable when J4 has finite mag- 
nitude. In the case of the ideal Kagome lattice, all the dimer 
coverings are degenerate up to the second order in A and the 
degeneracy can only be lifted by the third order perturbation 
processes as emphasized in Ref,i9. This is because the num- 
bers of the filled and empty triangles on the ideal Kagome 
lattice are fixed for any dimer covering due to the spatial ge- 
ometry of the ideal Kagome lattice (See Refi). This is, in 
fact, the reason why the 36-site-unit-cell VBS is the most sta- 
ble VBS state on the ideal Kagome lattice. Therefore, we may 
expect that the change in magnitude of J4 would affect the rel- 
ative importance between the second order and the third order 
processes. In the next section, we investigate the effect of fi- 
nite J4 on the ground state properties of the deformed Kagome 
lattice. 



III. VALENCE BOND SOLID PHASES OF THE 
DEFORMED KAGOME ANTIFERROMAGNET : THE 
EFFECT OF FINITE J4 

A. Dimer Series Expansion Study 

In order to refine the ground state energy estimation beyond 
the third order perturbation theory and confirm the results of 
the previous lowest order expansions, we develop a systematic 
scheme to compute the perturbation series using the linked 
cluster expansion method. In this approach, the ground state 
energy can be expressed as the sum of the contributions from 
a sequence of finite clusters-^-^^'. To carry out the series expan- 
sion, all the occupied links that make up the dimer covering 
are given an interaction strength J and all the other unoccu- 
pied links are given a strength A J. We apply the series ex- 
pansion method to compute the energy of the two competing 
VBS phases on the J4 -depleted lattice, J4 — 0, and the ideal 
Kagome lattice limit, J4/J ^ 1, up to the fifth order in A. 

The ground state energy per site of an infinite lattice can 
be written as a sum over cluster terms of the form, Eq /N = 
'^gc{g)£{g)- Here c{g) is the so-called lattice constant for 
the cluster g, which is the number of embeddings of the 
cluster g per lattice site. e{g) is the reduced energy of the 
cluster g, which is defined recursively as e(.g) — E{g) — 
J2g' c{g' / g)£{g')- Here E{g) is the ground state energy of the 
cluster g and the sum is over all subclusters g' of g. Since the 
reduced energy e{g) is equal to the sum of the contributions 
from all connected diagrams (i.e. the patterns of connected 
dimers) that span the cluster g, such diagrams only contribute 
beyond a certain minimum order in A, proportional to the size 
of the cluster. Thus, to obtain a perturbation series that is ex- 
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FIG. 4: (Color online) Topologies of the graphs (diagrams) relevant 
to the VBS-12 phase on the J4-depleted Kagome lattice that con- 
tribute to the ground state energy up to the fifth order in A. The thick 
dotted (pink) lines are the weak links on the hexagons while the thin 
dotted (black) links represent the weak links on the empty triangles. 



act up to a given order, it is sufficient to sum up the reduced 
energies of all the clusters up to the corresponding size. This 
gives us a systematic procedure for carrying out the expansion 
to higher orders. 

Another appealing aspect of applying the linked cluster ex- 
pansion method is that we can considerably reduce the com- 
plexity of the calculation by using the freedom in choosing 
clustersi212^ In a recent work of the dimer series expansion 
study on the ideal Kagome lattice^, the unoccupied links on 
each triangle are grouped together, which results in signifi- 
cant simplification of the computation and rapid convergence 
of the series for the ground state energy. In this scheme, for 
any dimer covering on the ideal Kagome lattice, every empty 
wd filled triangles constitute the elementary clusters. All the 
other clusters can be constructed by putting the elementary 
clusters together appropriately. Based on this specially de- 
signed clusters, it was shown that the series for the ground 
state energy can be determined up to the fifth order in A us- 
ing only five different graphs (or diagrams).'' Summing up the 
terms in the series, the ground state energy per spin of the 
VBS-36 phase on the ideal Kagome lattice is estimated to be 
—0.432 J. For comparison, we compute the ground state en- 
ergy of the VBS-12 phase on the ideal Kagome lattice. The 
same graphs can be used for the calculation and the summa- 
tion of the terms in the series up to the fifth order in A leads to 
the energy of —0.424 J. Note that even though the VBS phase 
with a 36-site unit cell is the ground state of the ideal Kagome 
lattice, the energy difference between the two VBS phases is 
quite small. 

Applying similar strategies to the J4-depleted Kagome lat- 
tice, we compute the series for the ground state energy of the 
two competing VBS phases. We first consider the VBS-12 
phase on the J4-depleted Kagome lattice. Here we group 
the three weak links on an empty triangle together and dis- 
tinguish them from the other weak links on the hexagons. In 
this scheme, only five different graphs (shown in FiglH con- 
tribute to the ground state energy up to the fifth order in A. 
The smallest cluster in Figllfa) is made of a single dimer and 
has the reduced energy, eo,a ~ ^jJ ■ FigHtb) and (c) show 
two clusters whose contributions start from the second order 
in A. The reduced energies of the corresponding clusters, i.e.. 
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£2,b and e2.c respectively, can be obtained after subtracting the 
subcluster energy Sp.a appropriately and are found to be 
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Finally, the resonance processes coming from the remaining 
two graphs in FigHfd) and (e) can lower the ground state en- 
ergy through the fourth order processes. The corresponding 
reduced energies of the two clusters, i.e., £4 ^ and £4 e can be 
obtained by the same procedure used earlier and using the fact 
that the subclusters are made of the graphs in Fig|4](a), (b) and 
(c). Since the reduced energies of all the other clusters that are 
larger in size than those in FigHfd) and (e), begin with O(A^) 
terms, the ground state energy can be determined exactly at 
least up to the fifth order in A with those five graphs shown in 
FigS After we determine the lattice constants of all the clus- 
ters, the ground state energy for the infinite lattice is found to 
be 

— = -- - - — A^ - 0.0097656A^ - 0.0079753A^ 
J 8 32 256 

This leads to the estimation of the ground state energy, 
-0.498 J for the VBS-12 phase. 

Following similar procedures, the ground state energy of 
the VBS-36 phase is obtained as 







= ----A2- 

J 8 16 



256 



-A-* - 0.0107422A* - 0.0055141A^ 



This leads to the estimation of the ground state energy, 
-0.473 J for the VBS-36 phase, which is higher than that of 
the VBS-12 state. This is consistent with the previous results 
of the perturbation theory and the VBS-12 phase is the ground 
state of the J4-depleted Kagome lattice. 

In summary, the VBS-12 phase is lower in energy on the 
J4-depleted Kagome lattice ( J4 = 0) while the VBS-36 phase 
has lower energy on the ideal Kagome lattice ( J4 = 1). Based 
on this result, we expect that there should be an energy level 
crossing between these two limits as the magnitude of J4/ J 
is changed from zero to one. In order to obtain the detailed in- 
formation about the level crossing, we use the self-consistent 
bond operator mean field approach to the two phases for an 
arbitrary strength of J4/ J in the next section. 



B. Bond Operator Mean Field Theory 

We study the relative stability of the VBS-12 and VBS- 
36 phases using the bond operator mean field theor y^^i^^'^^ 
for the deformed Kagome lattice with an arbitrary strength 
of J4/J < 1. In this formulation, the dimer singlet degrees of 
freedom are used as natural building blocks and the quantum 
corrections coming from the triplet fluctuations can systemat- 
ically be investigated^^'^^. 

Let us consider the two S = ^ spins constituting a dimer 
singlet, S/j; and S^. The Hilbert space is spanned by four 
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FIG. 5: (Color online) Ground state energies of the VBS-12 and the 
VBS-36 phases as a function of J4/ J (J is set to one) obtained from 
the bond operator mean field theory. Note that there is a level cross- 
ing around J4 = 0.97 J as indicated by an arrow. 



states that can be taken as a singlet state, \s), and three triplet 
states, \tx), \ty) and jt^). Then, the singlet and triplet boson 
operators are introduced such that each of the above states can 
be created from the vacuum |0) as follows: 
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To eliminate unphysical states from the enlarged Hilbert 
space, the following constraint needs to be imposed on the 
bond-particle Hilbert space: 



1, 



where a = x,y, and z, and we adopt the summation con- 
vention for the repeated indices hereafter unless mentioned 
otherwise. 

Constrained by this equation, the exact expressions for the 
spin operators can be written in terms of the bond operators: 
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where Ea/B-y is the third-rank totally antisymmetric tensor with 

Exyz — 1- 

We consider the antiferromagnetic Heisenberg model 
Hamiltonian, H — Jij^i ' Sj, where the exchange cou- 

plings of the neighboring spins are taken to be Ji = J2 = 
J3 = J > J4 (See Fig.[T]i. Utilizing the bond operator repre- 
sentation of spin operators, the Hamiltonian can be rewritten 
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FIG. 6: (Color online) Evolution of the triplet dispersion spectra of 
the VBS-12 phase as a function of J4/J (J is set to one). Note that 
the lower flat mode is fourfold degenerate while the upper one is 
doubly degenerate when J4/J = 1. 



solely in terms of bond particle operators. The hard-core con- 
straint among bond particle operators is imposed via the La- 
grange multiplier method. The quartic interactions between 
the triplets in the resulting Hamiltonian are treated by the self- 
consistent Hartree-Fock approximation. 

The computational procedure for the VBS-36 phase is basi- 
cally in parallel to that of the recent work on the 36-site-unit- 
cell VBS ground state of the ideal Kagome lattice^'. Here we 
present some details of the bond operator theory applied to the 
VBS-12 phase on the deformed Kagome lattice. Since the six 
dimers within the unit cell of the VBS-12 phase are equivalent 
as is obvious from the six-fold rotational symmetry of the or- 
dering pattern in Fig[3ja), the singlet condensate density (si_„) 
and the chemical potential /ii.„ can be set to be (sj „) = s and 
A*i,n=/^ in our mean field theory. Here i denotes the position 
of the unit cell and n indicates the dimer index inside the unit 
cell. The hard-core constraint on the bond-particle operators 
is imposed by adding the following Lagrange multiplier term. 

The quartic interactions between the triplet particles are de- 
coupled using the mean field order parameters and Q-y, 

where = {tl„J},ma) and = (ti,„Qij^,„a). We con- 
sider four mean-field variables, P.y and Q^, labeled by 7 = 1 
or 2. Here Pi and Qi denote the diagonal and off-diagonal 
triplet correlations between the dimers residing in the same 
unit cell while the other two variables, P2 and Q2 represent 
the triplet correlations between neighboring dimers in differ- 
ent unit cells. The above four order parameters Pi, Qi, P2, 
and Q2 together with s and /i are determined self-consistently 
by solving the coupled saddle point equationsi2i2^ 

In Fig. |5] we plot the ground state energies of the VBS- 
12 and VBS-36 phases as a function of < J4/J < 1. 
When J4 = 0, the ground state energy of the VBS-12 (VBS- 
36) from the bond operator mean field theory is — 0.496J 
(—0.470 J) which is close to the results from the series expan- 
sion study, —0.498 J (—0.473 J). As one increases the mag- 
nitude of J4, the ground state energy of the VBS-12 grows 



more rapidly than that of the VBS-36 phase and finally an en- 
ergy level crossing occurs around J4 = 0.97J. Beyond this 
point the VBS-36 phase becomes more stable and remains as 
the ground state up to the ideal Kagome lattice limit. 

The evolution of the triplet excitation spectrum for the 
VBS-12 phase shows interesting behavior as the magnitude 
of J4 changes. In Fig.|6] we plot the triplet dispersion spectra 
when J4 = 0,0.5 and 1, respectively. Note that the triplet 
dispersion of the VBS-12 phase on the ideal Kagome lat- 
tice ( J4 = J) shows completely flat spectrum. Flat disper- 
sion in the momentum space indicates the existence of local- 
ized eigenstates.^' In this case, the three neighboring dimers 
around an empty triangle support the localized eigenmodes. 
When J4 = J, the six dimers within a unit cell are all 'or- 
thogonal' to each other in the sense that each spin of a dimer 
is simultaneously connected to the two spins of its neighbor- 
ing dimers as shown in Fig. |3la). In this situation, the triplet 
particles can fluctuate only through the dimers of its neighbor- 
ing unit cell connected via empty triangles, and hence they are 
confined within that space. The three localized triplet modes 
coming from the three dimers around an empty triangle con- 
sist of doubly degenerate E modes and a nondegenerate A 
mode (in the standard convention of irreducible representa- 
tions in group theory.^-). Since the unit cell consists of six 
dimers, each E mode and A mode is doubly degenerate. In 
Fig. |6jc), the lower flat mode has fourfold degeneracy while 
the upper one is doubly degenerate. 

The fact that all the triplet dispersions are completely flat 
is an interesting characteristic of the triplet spectrum of the 
VBS-12 phase on the ideal Kagome lattice (note, however, 
that the true ground state is the VBS-36 phase on the ideal 
Kagome lattice). According to a recent work^', the triplet 
spectra of the VBS-36 phase on the ideal Kagome lattice con- 
tain some dispersive modes even though the number of flat 
modes is quite large there as well. In particular, the lowest 
flat band is degenerate with the second lowest dispersive band 
only at the zone center. As discussed in detail in ReffHHH 
the degeneracy at a particular momentum point reflects the ex- 
istence of noncontractible loop states and is topologically pro- 
tected. Therefore, the completely flat spectrum of the VBS-12 
phase on the ideal Kagome lattice means that all the flat modes 
are topologically trivial. This results from the fact that the size 
of every localized mode is not larger than the size of the unit 
cell and each localized eigenmode around an empty triangle is 
linearly independent. The condition for the emergence of the 
topologically nontrivial structure in dispersions is discussed 
in detail in ReflM (See also Refill) 

On the other hand, when the magnitude of J4 is reduced, 
the degenerate bands split and triplet particles spread over the 
whole lattice. In the case of J4 = 0.5 J that may directly be 
relevant to Rb2Cu3SnFi2, one can still see the gap between 
the lower four bands and the upper two bands. However, when 
J4 becomes vanishingly small, the gap structure disappears 
due to wide bandwidth of the triplet bands. 
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FIG. 7: (Color online) A local singlet resonance in the VBS-12 phase 
on the J4-depleted Kagome lattice. Here the dotted (red) square rep- 
resents the smallest flippable loop. 



IV. DISCUSSION 

In the case of the VBS phases with a big unit cell, a large 
number of degenerate dimer configurations may allow many 
singlet excitations below a spin gap. For example, the VBS-36 
order on the ideal Kagome lattice leads to the large degener- 
acy that comes from the combination of the broken transla- 
tional symmetry and the local resonance moves around flip- 
pable loops (See Ref 8 and Reffiol). Here the flippable loops 
mean a collection of alternating occupied and unoccupied 
links; when the roles of the occupied and unoccupied links 
on a flippable loop are interchanged, one hardcore dimer cov- 
ering can be continuously transformed to another one. The 
perfect hexagons and the pinwheel dimers form the smallest 
flippable loops in the case of the 36-site VBS phase on the 
ideal Kagome lattice.** The number of singlet excitations aris- 
ing from such considerations is still smaller than what was 
observed in the exact diagonalization studie&^i^^i^^i^ of small 
system sizes (up to 36 sites) for the spin- 1/2 antiferromag- 
netic Heisenberg model on the ideal Kagome lattice. It was, 
however, argued that the triplets can scatter to form additional 
singlet bound states in finite size systems and this may explain 
the exact diagonalization results'". 

On the other hand, the VBS-12 phase on the J4-depleted 
Kagome lattice does not break any translational symmetry and 
the local resonances from loop flips are the only sources of 
low energy singlet excitations. In Fig|7] we show the smallest 
flippable loop consisting of the four occupied and four unoc- 
cupied links, and its loop flip in the VBS-12 phase. Note that 
there is a change in the number of the empty and filled trian- 



gles after the loop flips. In fact, a single loop flip increases 
the number of the filled triangles by two. According to Eq.([T]| 
and Eq.(|2]i in the perturbation theory, this requires the energy 
cost of — 4Ae*^^) + 2Aei'^'' = ^J, which is much larger than 
the measured spin gap ~ 20K if the experimentally estimated 
J ~ 200K is used. Since the number of the filled triangles 
involved in a loop flip is proportional to the size of the flip- 
pable loop, the singlet excitations coming from a larger loop 
flip require even higher energy cost. Therefore, the VBS-12 
phase is quite stable against local singlet fluctuations and it 
seems difficult to produce "the singlet wave" whose existence 
is speculated in a recent experimental paper. We expect the 
lowest singlet excitation would be well above the spin gap, 
which can be contrasted with many singlet excitations below 
a spin gap in the case of the ideal Kagome lattice. 

In summary, slight lattice distortions can completely 
change the nature of the ground state of the quantum an- 
tiferromagnet on the Kagome lattice and this leads to the 
abrupt modification of the elementary excitation spectra in 
both the triplet and singlet channels. Even though Ji, J2, 
and J3 are not exactly the same in Rb2Cu3SnFi2, the mag- 
nitudes of these exchange couplings are very similar and sat- 
isfy Ji > J2 > J3 ^ J4. Interestingly, the VBS-12 phase 
studied in this work corresponds to the dimer covering that 
resides only on the Ji links. This suggests that the VBS-12 
phase on the Ji links is stable all the way to the limit where 
J2 and J3 become the same as Ji . Thus, we expect that our 
results are directly applicable to Rb2Cu3SnFi2. In particular, 
our predictions on the triplet and singlet spectra can be tested 
in the future neutron scattering and specific heat experiments 
on Rb2Cu3SnFi2. 
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